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In this paper we propose a new protocol to achieve coherent population transfer between two states 
in a three-level atom by using two ac fields. It is based on the physics of Stimulated Raman Adiabatic 
Passage (STIRAP), but it is implemented with the constraint of a reduced control, namely one of 
the fields cannot be switched off. A combination of frequency chirps is used with resonant fields, 
allowing to achieve approximate destructive interference, despite of the fact that an exact dark state 
does not exist. This new chirped STIRAP protocol is tailored for applications to artificial atoms, 
where architectures with several elementary units can be strongly coupled but where the possibility 
of switching on and off such couplings is often very limited. Demonstration of this protocol would 
be a benchmark for the implementation of a class of multilevel advanced control procedures for 
quantum computation and microwave quantum photonics in artificial atoms. 

PACS numbers: 03.67.Lx,85.25.-j, 03.65.Yz 


I. INTRODUCTION 

Preparation of a quantum system in a well defined 
state is an essential task in many branches of modern 
physics ranging from atomic and molecular physic^- 1 to 
quantum computation 2 . Techniques for transferring pop¬ 
ulation from a ground state |0) to a state |1) employ ei¬ 
ther Rabi cycling or adiabatic passage (APpl. Amongst 
these latter STIRAP is a three-level atom scheme where 
selective and faithful population transfer is achieved by 
ope rati ng with two resonant driving fields in A configura¬ 
tion. The advantage over Rabi cycling is the dramatic 
reduction of sensitivity to fluctuations of the parameters, 
at the expenses of a longer duration of the adiabatic pro¬ 
tocol. In more complex architectures semiclassical driv¬ 
ing fields are substituted by harmonic modes of a strongly 
coupled cavity, and tasks as preparation of photo ns w ith 
controlled amplitude, frequency and polarization 5 6 can 
be performed by vacuum-stimulated Raman AP (vSTI- 
RAP). 

In the last few years multilevel coherence in solid-state 
systems, from mesoscopic devices^ to atomic-like impu¬ 
rity stated, has been a fertile subject of investigation. In 
particular superconducting-based “artificial atoms’®^ 
are very promising since coherent phenomena proper of 
the microscopic realm have been demonstrated on the 
mesoscopic scale. With respect to their natural counter¬ 
part, artificial atoms offer the advantage that composite 
structures can be fabricated on a single chijP-- : , which al¬ 
lows manipulation of single photons at GHz frequencies 
opening the new scenario of microwave quantum photon- 
icsP. Moreover new architectures could be implemented 
with both larger coupling:!^ and a larger degree of inte¬ 
gration than their atomic counterparts. 

In the last few years several theoretical proposals-®^ 
and experiments^] h av e dealt with multilevel coher¬ 
ence in artificial atoms. Distinctive features of s uch sys- 
tems are the effectiveness of parity selection rules 16 18 21; 


which together with the presence of strong 1// nois d 26 * 2 ^ , 
impose constraints on the available control. Therefore 
new protocols for manipulating the coherent dynamics 
must be tailored for such systems. Their design requires 
that large couplings allowing for efficient control are com¬ 
bined with protection from nois^Il. 

In this paper we present a new protocol to achieve 
coherent population transfer between the two lowest ex¬ 
cited states of a three-level atom by using two ac fields. 
The key difference with standard STIRAP, where ac 
fields must be switched on and off in a counterintuitive 
sequenc^, is that one of the fields is kept always-on, its 
amplitude being constant during the protocol. Opera¬ 
tions require phase modulation, and for this reason we 
call the protocol cSTIRAP (chirped STIRAP). Sweeping 
the frequency of a single classical driving field across the 
resonance is a standard way to switch on and off Rabi 
oscillations, thereby one may think to rephrase STIRAP 
accordingly, but this is not the case. Indeed coherent 
population requires destructive interference of the two 
Heldd^. This is guaranteed by cSTIRAP, which thereby 
solves a non-trivial control problem, its experimental 
demonstration in artificial atoms being by itself an im¬ 
portant proof of principle of advanced three-level control. 
Even more interestingly, cSTIRAP could apply to archi¬ 
tectures where “artificial atoms” are coupled to quan¬ 
tized modes, electromagnetic or nanomechanical, where 
strong coupling is achieved by non-switchable hardware 
elements keeping the interaction always-on. The proto¬ 
col we propose possesses certain advantageous distinctive 
characteristics: (i) it works with reduced available con¬ 
trol, as always-on held, (ii) it operates with nearly reso¬ 
nant fields, reducing the operation time; (iii) it may rely 
on better techniques to control the phase of microwave 
circuits, (iv) it is cyclic. 

The paper is organized as follows. In Sec|ll] we in¬ 
troduce the model Hamiltonian and briefly review stan¬ 
dard implementations of coherent population transfer in 
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two and three-level atoms. In Sec. uni we illustrate the 
new protocol discussing in Sec.[IV]the robustness against 
parametric fluctuations and in Sec. |V| decoherence effects. 
Finally, in Sec |VI[ along with the conclusions, we will dis¬ 
cuss the comparison of cSTIRAP with other protocols for 
population transfer operated by frequency chirps. 


II. COHERENT POPULATION TRANSFER IN 
LAMBDA ATOMS 

In two-level systems coherent population transfer 
|0) —>■ |1) by AP is performed by shining a direct coupling 
field whose detuning is swept through the resonance at 
the Bohr frequency of the transition. Common examples 
are Rapid AP (RAP) or Stark Chirped RAP (SCRAP)fS. 

In three-level systems population transfer may be 
achieved in absence of direct coupling, via a third linkage 
state 12), coupled to |0) and |1) by a pump field at fre¬ 
quency ojp cs E 2 — E 0 and a Stokes field at w s ~ E 2 — Ei , 
respectively. In particular the Lambda configuration de¬ 
picted in the top inset of Fig. [Tj will be considered in this 
work. Since |2) is usually short-lived, one of the goals 
of coherent techniques is to use |2) but never populate 
it. This is achieved in a very efficient and elegant way 
relying on destructive interference 28 . The Hamiltonian 
in the rotating wave approximation, in the basis of the 
bare states {|0), 11), 12)}, is expressed in a doubly rotating 
frame as 


H = 


0 

0 

2 Llp(t) 


0 

m 


Sp{t) 


(i) 


where flfe(t) with k = p,s are the Rabi frequencies of 
the pump and the Stokes fields, which are detuned by 
S p := E 2 — Eq — ajp and S s := E 2 — Ei — u> s respectively. 
A key quantity is the two-photon detuning, defined as 
6 :=5 P - S s . 

Conventional STIRAP relies on the fact that at two- 
photon resonance, S(t) = 0, an instantaneous eigenvector 
with zero eigenvalue eo = 0 exists given by 


= n s (t) | 0 ) - n p (t) | 1 ) 

yJW s (t)+W p (t) 


( 2 ) 


two-photon effective coupling |0) 0 |1). STIRAP has 
been observed in a variety of physical systems^. The 
two-photon character of population transfer, and the fact 
that the protocol is maximally efficient with fully reso¬ 
nant fields, S = S s = S p = 0 is the key for interesting 
applications with quantized fields. 

Another three-level technique, Raman Chirped Adi¬ 
abatic Passage (RCAP)pD, uses instead phase modula¬ 
tion. Population transfer is achieved by two far off- 
resonanc e ch irped laser pulse sweeping through resonance 
(see Sec. VI). Unlike conventional STIRAP, two-photon 
resonance is not kept during the whole process, causing 
a transient population of state |2) to appear. The lat¬ 
ter in principle can be made small by accurate tuning of 
parameters. 


III. COHERENT POPULATION TRANSFER 
WITH AN ALWAYS-ON FIELD 


In this section we will address the problem of achieving 
| 0 ) -A | 1 ) population transfer subject to two constraints, 
namely (a) keeping the population of | 2 ) small and (b) 
operating with a reduced control, in particular with one 
of the fields, for instance the Stokes one, kept always on, 
Ll s (t) =: O 0 7 ^ 0. Naively one could suppose that sweep¬ 
ing the detuning S s (t) could allow to effectively switch 
on and off fi s , allowing again for conventional STIRAP. 
However this is not the case because coherent population 
transfer requires that the two-photon resonance condi¬ 
tion, S = 0, is kept while sweeping 5 s (t), to ensure de¬ 
structive interference. In what follows we will seek for a 
solution allowing to achieve complete population transfer 
by properly shaping the detunings. 

First of all when one of the fields is always on, the 
Hamiltonian <[T|) for t -A ±00 is not diagonal in the bare 
state basis. In order to approximate asymptotically the 
desired target state 11 ), necessarily at the end of the pro¬ 
tocol we must have S s ^> LIq. If we take detunings shaped 
as shown in FigjT] which are given by 


S s (t) 

5 p {t) 


1 

2 


hs Uq 



KsS s (t) 


■ tanh 


t + T 
Tch 


( 3 ) 


It is called dark state since population is confined in the 
“trapped subspace” {|0), 11)}, despite of the fact that the 
two fields excite both the 0 —► 2 and the 1 -a 2 transi¬ 
tions. The key phenomenon preventing population of |2) 
is destructive interference between the amplitudes cor¬ 
responding to the two absorption patterns 1 -^. Conven¬ 
tional STIRAP consists in letting the dark state evolve 
adiabatically from \D(— 00 )) = |0) to |D(oo)) = |1). This 
is achieved by shining pulses flk(t) in a “counterintu¬ 
itive” sequence, the Stokes at first and then the pump. 
An important characteristic of STIRAP is the fact that 
AP is operated when both fields are on, determining a 


the desired asymptotics is ensured by hs 1, i.e. the 
protocol must start and end with “far detuned” lasers. 
The modulation ©> has the important characteristics 
that at least for part of the protocol S(t) = 0 (Fig. [l]). 
During this phase a Stokes-induced Autler-Townes (AT) 
splitting opens. Although an exact adiabatic dark state is 
not available for population transfer, we will argue later 
that keeping S ss 0 allows to minimize the transient pop¬ 
ulation of 12). 

The population transfer mechanism is better under¬ 
stood studying the evolution of the instantaneous eigen¬ 
values and eigenvectors of the “Stokes” Hamiltonian, ob- 
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FIG. 1. (color online) Main figure: single (coloured lines) and 
two-photon (dotted line) detunings in Do units. Top inset: 
Three level Lambda system. Bottom inset: Zoom of the sin¬ 
gle and two-plioton detunings (solid lines), plotted together 
with the Stokes eigenvalues (dashed lines) of Eq. © showing 
the appearance of a dynamical Stokes-induced AT, which is 
switched on and off by modulation of <5 S . 


tained setting to zero the pump field in Eq. 0 


H s {t) = 


0 0 0 

0 S(t) 

0 S p (t) 


(4) 


Here the Rabi frequency has been taken real with no loss 
of generality. The Stokes Hamiltonian acts non-trivially 
only on the {11), 12)} subspace, yielding the asymptotic 
states 


|s+(-°o)> - I 2 ) IM+°°)> - I 1 ) 

|s_(-oo)) ~ |1) |s-(+oo)) ~ |2) 


The ’’Stokes eigenvalues” display the presence of the AT 
splitting during the protocol (Fig. [lj bottom inset) 


s 0 — 0, 


s± = S + 


s s ±ysf +^ 0 

2 


( 6 ) 


During this AT phase S s is swept across the resonance 
swapping |1) o |2). 

Using detunings Eq.(0 with Kg > 1 the pattern of split 
instantaneous eigenvalues s±(t) is crossed twice by the 
eigenvalue So = 0, as shown in Fig. 0a). Crossings occur 
at times ±t c when s±(t) = 0, i.e. 4S(t c )Sp(t c ) = Ug. 
In these conditions the system prepared in \ip(— oo)) = 
|0) remains of course in this state, passing through the 
crossing. Population transfer is achieved by applying a 
pump pulse with finite area reaching its peak value close 
to the second crossing, t = t c . For instance, employing a 
a Gaussian pulse, we have 


£l p (t) = Ktt 0 e v T ) 


(7) 


The behavior is understood in terms of the instanta¬ 
neous eigenenergies of the full Hamiltonian Eq. 0. In 
particular the pump pulse lifts the degeneracy between 
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FIG. 2. (color online) (a) Eigenvalues of the Stokes Hamil¬ 
tonian of Eq. (pi) ) (dashed lines) and of the full Hamilto¬ 
nian of Eq. (solid lines) in Do units. The red thick 
line is the instantaneous energy of the system adiabatically 
driven from |0) to |s+) — |1) through the opening of the 
avoided crossing generated by the pump pulse (dotted line) 
at time t = t c . (b) Population histories (red, blue and green 
lines) from the numerical solution of the Schrodinger equa¬ 
tion, for f l s (t) = Uo,UoT = 40, hs = 10, ks = 1.2 and k = 1, 
Tch = 0.6 T, showing complete population transfer by cSTI- 
RAP. For these parameters the adiabatic approximation (gray 
curves) fully agrees with the exact solution. Inset: the exact 
population P 2 of the excited state (green solid line) is small 
at any time of the protocol, as can be estimated by Eq.(|8| 
(thin line). The dashed line refers to the approximation of 
Eq. (|A4l). 


So and s + turning their crossing into an avoided cross¬ 
ing [Fig0a)]. The adiabatic connection corresponding 
to s+ yields eventually the desired population transfer, 
|0> -A |s+(+oo)) ~ |1). 

We remark that population transfer depends only on 
the presence of a crossing between Stokes eigenenergies 
s + and So and on the fact that the adiabatic approxima¬ 
tion is valid. In this regime the precise shape of the pulses 
is not relevant. Therefore the protocol is robust against 
imperfections in the control. From the physical point of 
view it is worth stressing that the pump pulse triggers 
AP by a two-photon process. The distinctive feature of 
our proposal is that this two-photon effective coupling 
is obtained with both quasi resonant pump and Stokes 
fields. This ensures large efficiency for rather small pulse 
duration. We mention that during its switching on Q p 
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could in principle trigger unwanted transitions |0) —> |2), 
which are however suppressed by the Stokes-induced AT 
splitting and the two-photon resonance condition. A sim¬ 
ilar phenomenon occurs in standard STIRAP, where it is 
called the Stokes-induced EIT (electromagnetically in¬ 
duced transparency) phas^. 

Summing up cSTIRAP can be described in the lan¬ 
guage of Ref. 1 as a five stages protocol, with suc¬ 
cessive far-detuned, Stokes-induced AT, Stokes-induced 
EIT, two-photon AP and again far-detuned phases. In 
what follows we will see that the other important require¬ 
ment, namely that population of |2) should be minimal 
at all times, is also fulfilled. This requirement is neces¬ 
sary in order to prevent unwanted decay processes likely 
to occur in real physical systems, where |2) is often un¬ 
stable. 

We estimate P 2 = |(2|7)| 2 by adiabatic elimination. 
The standard procedure formulated in the bare basis 
yields the state | i^ae) = c o|0) + c i|l) ( see App. 0 
First order corrections yield a leakage from the subspace 

{|0),|1>} given bjEU 


P 2 W 


UpCo —f- floci 

25 P 


2 


( 8 ) 


that can be made very small, as shown in Fig|2])b) , 
which also shows that this approximation works very 
well. A better approximation is obtained by working in 
the Stokes basis [see App. a and Fig. [2][b)}, but Eq.([8]) 
has a simpler analytic form, allowing to write a figure 
of merit for the parametric dependence of leakage during 
the protocol. A simple choice is to consider leakage at 
the crossing s + = 0 
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FIG. 3. (color online) Upper panels: STIRAP efficiency vs 
the relative peak amplitudes of the fields (left panel, where 
Ks = 1-2) and to the relative detunings (right panel, where 
k = 1), for various degrees of adiabaticity (curves: floT = 40 
(red), 20 (blue), 10 (brown) from higher to lower efficiency). 
For /tf > 1, and provided adiabaticity is good, the system has 
a very slight sensitivity to these parameters. Lower panels: 
sensitivity of the efficiency to unwanted transient population 
of 12) accounted for by a finite lifetime T 2 =: I/T 2 (cf. Eq[K)]). 
The insets of panels (a) and (c) are zooms of the corresponding 
main figures showing how, in the presence of a non-vanishing 
r 2 , the efficiency improves with increasing ft. 


eigenvalues do not cross, and adiabatic dynamics leads 
to a final population entirely in |0). 


IV. SENSITIVITY TO PARAMETERS 


P2{t c ) - = KS . 1 f( K ) (9) 

Op 

Here /(ft) is a monotonically decreasing function of the 
ratio of the Rabi peak amplitudes ft. This qualitative 
behaviour is confirmed by numerical simulations shown 
in Fig. [3j where the efficiency is plotted versus relative 
magnitude of the amplitudes (ft, left panel) and of the de¬ 
tunings (ft^, right panel), both in the absence (top panel) 
and the in presence (bottom panel) of a finite lifetime 
T2 = T/2 of the intermediate state |2) (see section |Tv| for 
a model). It is seen that efficiency increases with increas¬ 
ing ft as an effect of a larger avoided crossing at s + = 0. 
Moreover increasing ft reduces the transient population 
of 12), as given by the figure of merit Eq.([9]). This is 
seen by comparing the insets of the left panels of Fig. [3] 
the positive slope of the sensitivity in the presence of a 
finite T2 [Fig. [3jc)] cannot be explained as an improve¬ 
ment in adiabaticity, since this slope is not present in the 
ideal case [Fig. |3](a)]. Therefore, it can only be caused 
by a reduction of P^- Population transfer occurs only 
for Kg > 1 as shown in Fig. |3][b),(d). In particular, for 
Kg = 1 we have S(t) = 0 and Eq.Q applies, showing that 
an always on Stokes field would produce a return of the 
population to the initial state. For Kg < 1 the Stokes 


The efficiency of cSTIRAP is not very sensitive to 
slight deviations of relative amplitudes ft and detunings 
Kg of the pulses, provided adiabaticity is kept. This is 
shown in Fig. [3j results in the lower panels allowing to 
fix convenient values k = 1, Kg = 1.2 and f 1 0 T = 40, 
which we will use hereafter. 

As in conventional STIRAF^, the most critical feature 
is the parametric sensitivity to stray detunings. Here we 
discuss this issue, which is als o resp onsible for decoher¬ 
ence due to low-frequency noix 21 

The physics is understood recalling the picture of con¬ 
ventional STIRAP, where two kind of errors ernergd®. 
“Bad projection” errors, due a bad choice of the 
pulse shape and timing, may lead to the wrong target 
state. “Bad adiabaticity” errors induce leakage from the 
trapped subspace, nonadiabatic transitions surely occur¬ 
ring when the so called “global condition” 1 

is not met. Both kinds of errors are also triggered by 
fluctuations induced by an environment (see SecjV|. For 
cSTIRAP we verified that large enough again guar¬ 
antees adiabaticity (Fig. [ 3 ]). In this regime a strong asset 
of cSTIRAP is that it is not affected by bad projection 
errors in the far-detuned phases, since final eigenstates 
in the rotating frame are nondegenerate. 
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However since the efficency of cSTIRAP depends on 
the structure of crossings of the eigenvalues of the Stokes 
Hamiltonian, it may be affected by stray detunings dur¬ 
ing the protocol. A further drawback comes from the fact 
that the state carrying population in cSTIRAP, while 
taking advantage from destructive interference, it is not 
an exact dark state as in Eq.([2]), since the condition 
8(t) = 0 does not hold true. This is a potentially im¬ 
portant source of error for cSTIRAP since it also may 
determine a nonvanising population of \2) at intermediate 
times. Sensitivity to detunings is conveniently studied by 
the non-Hermitian Hamiltonian: 


H(t\{8 k }) ^ H(t\{8 k }) +iT 2 \2)(2\ 


( 10 ) 


Using a sufficiently large T 2 > 1/T guarantees that tran¬ 
sient population of |2) decays elsewhere (e.g. in a con¬ 
tinuum), yielding a lack of normalization at the end of 
the protocol. Therefore the resulting efficiency Pi(tf) 
is a figure of merit embedding the requirement that \2) 
should be never populated. 

The Hamiltonian (10), where only the dependence on 
detunings is emphasized, accounts for the effect of stray 
components by letting 


Sk(t) -» 8 k (t) + S k , k = s,p 

S(t) S(t) + ~S, 5 := 8 P — S s ' ' 

Stray detunings may describe very slow phase fluctua¬ 
tions (at frequencies -C 1/T) of the driving fields. Phys¬ 
ically in solid-state devices they describe energy fluctua¬ 
tions due to coupling to an environment (see Sec fV| a nd 
Ref. ETT) whose power spectrum has 1 //“ behavior^. In 
what follows we describe the detrimental effects they pro¬ 
duce and the limitations they determine. 

The efficiency of the protocol versus stray detunings is 
shown in Fig. [ 4 ] The colour map shows T\ (t) for P 2 = 0 
at the end of the protocol, t = tf. Lines refer to finite 
T 2 = 1/T, which determines a reduced value of P\{tf) 
since a nonvanishing population P 2 (t) would decay out¬ 
side the system. It is seen that the efficiency is large 
in a whole region around the center of the plot (absence 
of fluctuations, 8 S = 8 P = 8 = 0), showing the stability 
of the protocol. The failure of cSTIRAP in the region 
of larger detunings is analyzed in the App. [C] Here we 
mention that in the first and in the third quadrants of 
Fig. a failure is due to “bad projection” errors, i.e. the 
system may evolve along an adiabatic linkage leading to a 
wrong target state. Instead deep in the second quadrant 
the problem is “bad adiabaticity” due to an insufficient 
pump-induced two-photon avoided crossing. 

Concerning sensitivity to r, notice that the convenient 
delay is implicitly set by the choice of fl p (t) being maxi¬ 
mal at the second crossing time, Eq. |7]). We have checked 
that in these conditions the protocol is stable against 
deviations from the delay and the detailed pulse shape 
used in this work, provided they are not too large. More¬ 
over it is worth stressing that the protocol we propose in 
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FIG. 4. (color online) The color map describes the efficiency 
of ideal cSTIRAP, with lb = 0 vs fluctuations of the de¬ 
tunings. In the brightest area we have Pi(tf) > 0.9. Lines 
refer to lb = 1/T and delimit the Pi(tf) > 0.9 (most inner 
region) and the Pi(tf) > 0.8 areas. We use the same pa¬ 
rameters as in Fig.j2jb) , which guarantee that in absence of 
fluctuations, <5 S = 8 P = 5 = 0, adiabaticity of the protocol is 
strong. The extension of the regions of large efficiency deter¬ 
mines the single-photon linewidths (in this case A<5 S ) and the 
two-photon linewidth A<5. 



the “ideal” detunings case, while being physically satis¬ 
factory, is not an optimal solution in the mathematical 
sense. Therefore we expect further improvement by tack¬ 
ling the problem with Optimal Control Theory. 


V. DECOHERENCE 


A further important source of errors in STIRAP is de- 
coherenc^H, especially in solid-state artificial atoms. We 
discuss some qualitative aspect in this section. A key 
asset of conventional STIRAP is that while spontaneous 
decay from \2) may be large (decay time larger than the 
duration of the protocol), the phenomenon is supposed 
to have small impact as long as \2) is depopulated. This 
holds true also for cSTIRAP, as seen from the results 
for P 2 7 ^ 0 presented in the last Section. Markovian de¬ 
phasing in STIRAP has been studied in detaiP^ and its 
detrimental effect, namely leakage from the trapped sub¬ 
space due to the weakening of destructive interference 
phenomenon, has been elucidated. It has been shown 
that strong Markovian dephasing is tolerated, as long 
as it does not affect the two levels of the trapped sub¬ 
space. More complete studies of the effects of quantum 
noise in driven systems have pointed out that in solid 
state implementations of three-level artificial atoms the 
main effect is due to decay processes within the trapped 
subspac4-0. Other decoherence channels emerging in the 
Born-Markov approximation, namely the relation of rates 
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to the detailed spectral density of the environmental and 
the possible drive-induced absorptiorPS are less relevant. 

On the other hand, unlike their natural counterpart, 
artificial atoms implemented by solid-state nanodevices 
suffer from low-frequency noise 2 '. This drawback may be 
compensated by the ease of producing large couplings on 
chip, the tradeoff between protection and addressability 
being the central design issue. The effect of low-frequency 
noise in STIRAP has been discussed in Ref. [20,; where 
its interplay with Markovian n oise and the role of de¬ 
vice design were also addressecP 1 ^!. The extension of 
this detailed analysis to cSTIRAP is beyond the scope of 
this paper, but general features pointed out in the above 
works together with the results of the last section, allow 
to draw a physical picture which can be used as a guide 
for device design. 

We assume that low efficiency may be determined by 
by decoherence leading to detrapping from the {|0), 11)} 
subspace and by failures of the adiabatic approxima¬ 
tion also leading to unwanted population of |2). The 
simplest model encompassing these main features is to 
account for decay of |2) in a continuum due to quan¬ 
tum noise (T 2 ) and to account for dephasing as due 
to low-frequency (classical) fluctuations of relevant pa¬ 
rameters. That is we consider the Hamiltonian Eq. 0 
sup plem ented by the non-Hermitian term appearing in 
Eq.(10 ), and let Skit) —> Skit) + Skit), for k = s,p, and 
f2fc(t) —» ilk {t) + ilk{t), where Skit) and ilk it) are clas¬ 
sical stochastic processes. In artificial atoms such fluc¬ 
tuations stem physically from noisy external bias fields, 
which induce fluctuations of the energy splittings of the 
device (determining Sk’s) and of the operator coupling 
to the held (yielding fife’s). The efficiency is obtained 
by averaging over such fluctuations -Pi(t|{tffc}, {f2fc}), at 
the end of the protocol. In cases of interest, as for 1/f 
noise, the average can be est imate d in the quasistatic 
(or static-path) approximatiorP^ 2 3 It amounts to sub¬ 
stitute stochastic processes by random variables with a 
suitable Gaussian distribution, which physically accounts 
for sample to sample fluctuations of parameters. Results 
of the last section suggest that stray HkS hardly affect 
the efficiency, whereas the effect of the distribution of 
Sk’s can be important. This effect can be read off in 
Fig|4j which shows that for reasonably small fluctuations 
there is a region where still large efficiencies are allowed. 
Successful cSTIRAP requires that fluctuations of energy 
levels (i.e. detunings) are smaller than the linewidths. In 
analogy with the analysis of Ref 121] we expect that the 
condition of large efficiency depends on the bandstruc- 
ture of the device at the bias point. Indeed depending on 
the device and on the noise source, fluctuations of the two 
splittings (detunings) are either correlated or anticorre- 
latecP^, namely they are described by lines with positive 
or negative slope in Fig. [dj A figure of merit is the ra¬ 
tio Si /as between the two-photon linewidth of STIRAP, 
corresponding to the width of the large efficiency region 
in the proper direction in Fig. [4j and the variance as of 
the fluctuations of the two-photon detuning. 


VI. CONCLUSIONS 

In this paper we have proposed a new protocol 
which extends conventional STIRAP. Coherent popula¬ 
tion transfer is achieved with reduced available control, 
namely one of the field is kept always on. This pro¬ 
cedure is suited for applications in artificial atoms and 
can be advantageous in integrated atom-cavity systems 
architectures, where couplings to quantized modes are 
implemented by non-switchable hardware 12 ^, and may be 
manipulated in this way for applications to microwave 
quantum photonics 13 . In this respect it may be useful 
that cSTIRAP can be repeated cyclically since popula¬ 
tion histories are invariant when Sk —> —Sk, allowing the 
protocol to work as well in the reverted detunings con¬ 
figuration. 

The protocol leverages on the fact that in the mi¬ 
crowave realm external fields have a phase which can 
be usually controlled better than for sources at optical 
frequencies. In particular frequency can be modulated 
more accurately allowing direct time-dependent control 
of the detunings, instead of the induced Stark shifts used 
in genuine atomic system^ 1 -. Moreover in solid-state arti¬ 
ficial atoms, e.g. superconductor based, detunings can be 
independently modulated by external voltages and fluxes. 

Manipulation of detunings is the basis of other coher¬ 
ent transfer protocols like RCAF^. The essential dif¬ 
ference between standard RCAP and cSTIRAP is that, 
owing to the fact that the Stokes field is always-on, our 
protocol involves a dressed state in the AP phase (see 
Sec. Ill), whereas in the former AP occurs between bare 


states. Therefore while in RCAP the avoided crossing 
is due to the two-photon coupling of two far detuned 
dispersively coupled fields, in cSTIRAP AP takes place 
via destructively interfering resonant fields. This renders 
more robust the protocol, which achieves large efficiency 
for rather small pulse duration. On the other hand the 
analogy with RCAP, as well as the discussion of Sec. |Vj 
suggests that also cSTIRAP may be resilient to phase 
noise and to low-frequency noise in nanocircuits offering 
adva ntag es in quantum state processing with artificial 
atoms—. 

STIRAP is also the basis of other protocols as prepa¬ 
ration of superposition^ transfer of wavepacketpD, ma¬ 
nipulation of photons and quantum gates^, with still 
unexplored potentialities for quantum information and 
quantum control. Therefore demonstration of cSTIRAP 
is a benchmark for a class of multilevel advanced control 
protocols in artificial atoms. 


Appendix A: Adiabatic elimination of state |2) 

In order to estimate the population of |2) we start 
from the usual adiabatic elimination in the bare ba¬ 
sis. The Schrodinger equation id t \tp) = H\ip), with the 
Hamiltonian Eq.Q, is written for the components of 
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FIG. 5. (color online) Population histories for =: 

flo,^oT = 40, hs = 10, Ks = 1.2, k = 1 and r c h = 0.6 T. 


c_ ~ 0, we obtain c_ = — (fi_/2s_)co, where Cl± = 
f2 p [l + 4(<5 — Sqz/Ilo) 2 ]” 17,2 . Substituting the Ansatz for 
c_ into the Schrodinger equation yields an effective 2x2 
Hamiltonian, which in the {|0), |s+)} basis reads: 


H 2s = 


fP n + ~ 



(A3) 


which yields the leakage to |2) in the form 


P2 


- i +|2 

-—c 0 o 2 +c + aj | 
2s- 


(A4) 


As it is seen from Fig. [2jb) (dashed line) the result re¬ 
produces the numerical solution, but it does not yield a 
figure of merit as simple as Eq. [8| 


\4>) := Y^i=o c il*)- Assuming c 2 — 0 one finds 


Appendix B: Always-on pump field 


c 2 = - 


f2pCo —f- HgCr 

2 T P 


(Al) 


This expression of c 2 is substituted in the Schrodinger 
equation yielding a two-state problem governed by the 
effective Hamiltonian 


H 2 (t) = 


n s g p 

4 S p 4<5 p 

Q s £ip c 

4 S p 4S p 


(A2) 


Now assuming the validity of the adiabatic approxima¬ 
tion, Co and ci are approximately given by the instanta¬ 
neous eigenvectors of H 2 (t). In particular we consider the 
state corresponding to the preparatio n | il> (tj)) = |0), and 
we can estimate P 2 = \c 2 \ 2 from Eq. (Al). The analytic 
result is shown in Fig. [2jb) , thin solid line in the inset, 
and it yields good agreement with the numerical curve. 
The analytic expression, though easy attainable, is cum¬ 
bersome. Insight in the parametric dependence can be 
gained by evaluating leakage at t = t c : 


PHI- Ks ~ 1 («~v'k 2 +4) 2 

2 cj ~ 4 +(k + 7)?T4) 2 

which is Eq. Q. We remind that adiabatic elimination 
yields coarse grained amplitudes and it is a priori en¬ 
forced by large single-photon detunings. Remarkably the 
result obtained from Eq.( |Al[ ) is accurate for the whole 
procedure, even if in part of the protocol the condition 
S p flk is not met. This is due to the fact that the pop¬ 
ulation of 12) is always small, either because the regime 
is dispersive or because there is destructive interference. 

Corrections in the regime where S p (t ) < can 

be fully taken into account if adiabatic elimination 
is carried in the representation of the Stokes eigen¬ 
states. We write the Hamiltonian Eq.Q in the basis 
{|0>, |s+), ]«_>}, given by |s±) = af |1) + af\2). By ex¬ 
pressing j'i/j) = cq|0) + c+|s+) + c_|s_) and assuming 


We can seek for a protocol dual to the always-on Stokes 
field by making the following substitutions, t c —► —t c , 
S P W 8 S , Sip ^ fl s . The population histories are shown 
in Fig. [5] and differ somehow from those of Sec [Hi The 
point is that the system is prepared in 10), which in this 
case is not an exact eigenstate of the initial Hamiltonian. 
As a consequence Rabi oscillations of small amplitude 
appear in both Pq and P 2 . They can be substantially 
reduced by increasing the initial value of the pump de¬ 
tuning. Stray population may appear in the intermediate 
state 12) also due to adiabatic population transfer at the 
avoided crossing, and can be minimized by adjusting pa¬ 
rameters as suggested by Eq. (8j 


Appendix C: Failure of STIRAP at large detunings 

We now analyze the dynamics in the regions of Fig. [4] 
where cSTIRAP fails. As mentioned in Sec. HVl when en¬ 
ergy levels have infinite lifetime, failure of the protocol is 
due to two kind of errors, namely ” bad adiabaticity” and 
’’bad projection” 3 . While in the former case, the protocol 
fails because the avoided crossing produced by the fields 
is insufficient to guarantee adiabaticity, in the latter case 
the system is projected onto the wrong eigenstate of the 
Hamiltonian. Errors mainly occur during the AP near 
the point at t = t c where Stokes eigenstates cross. An 
efficient protocol requires for the probabilities of Landau- 
Zener transitions between such states that 7o-> s _ 1 

and 1 — 7o->s + 1, which is not always met for finite 

stray detunings. 

A qualitative picture of how cSTIRAP possibly fails 
due to stray detunings is offered by the patterns of the 
instantaneous eigenvalues of the full and of the Stokes 
Hamiltonians, in the darker regions of the three (a-c) 
quadrants of Fig. [4] Examples of these patterns are plot¬ 
ted in Fig. [fija-c). 























FIG. 6. (color online) (a-c) Instantaneous eigenvalues of the 
Stokes (dashed lines) and complete (solid lines) Hamiltonians, 
in the dark regions of quadrants (a-c) of Fig. [4] in units of 
H 0 . (d) Population histories corresponding to quadrant (d) 
of Fig. [4] for F2 = 0, showing that, even if Pi{tf) is nearly 
one, the protocol suffers of large transient population of |2) 
(green line). In each panel the label ( 6 /Q 0 , <5p/Ho) indicates 
the value of the stray detunings. 


sequent swap |s_) —» |2) leads to a wrong target state. 
Deep in quadrant (b), the protocol suffers from a sort of 
final “bad projection” error: the second crossing occurs 
at negative times , where Cl p ~ 0 and the correspondent 
transition becomes diabatic. This yields | ip(t)) ~ |0) at 
all timeJ®. Deep in quadrant (c) cSTIRAP fails when 
the configuration of detunings renders the pump-induced 
avoided crossing insufficient. In this case the problem is 
“bad adiabaticity”, Zener tunneling inducing unwanted 
transitions to the state adiabatically evolving towards |0). 

Finally, deep in the quadrant (d) the configuration of 
detunings is such that the two “mixing” phases of the 
protocol are inverted. Indeed the Stokes-induced AT 
splitting becomes relevant only after the second cross¬ 
ing, which in the ideal case would have produced the 
two-photon AP. Therefore fl p partially injects popula¬ 
tion into |2). At later times, in the Stokes-AT phase, 
this population is swapped to |1). Although the final 
state is correct (cf. the large efficiency in Fig. [4]), in the 
presence of decay F 2 7^ 0, occupation of |2) at intermedi¬ 
ate times suppresses the efficiency [see Fig. [b](d) and the 
solid lines of Fig. [4]. 


In the region deep in quadrant (a) of Fig. [4] detunings 
are such that the first crossings of the Stokes eigenen- 
ergies occurs at positive times , i.e. when fl p is already 
on [Fig. [6ja)] . Therefore |0) and |s_) mix, originating 
a sort of initial “bad projection” error. Then the sub- 
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